Suppression of finite-size eff'ects in one-dimensional correlated systems 
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We investigate the effect of a non-uniform deformation applied to one-dimensional (ID) quantum 
systems, where the local energy scale is proportional to gj = [sin(j7r/A'')]™ determined by a positive 
integer m, site index 1 < j < A'^ — 1, and the system size A*'. This deformation introduces a smooth 
boundary to systems with open boundary conditions. When m > 2, the leading 1/N correction to 
the ground state energy per bond e^^' vanishes and one is left with a correction, the same 

as with periodic boundary conditions. In particular, when m — 2, the value of Eq^' obtained from 
the deformed open-boundary system coincides with the uniform system with periodic boundary 
conditions. We confirm the fact numerically for correlated systems, such as the extended Hubbard 
model, in addition to ID free-Fermion models. 
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I. INTRODUCTION 

The periodic boundary conditions (PBC) are often 
more convenient than the open boundary conditions 
(OBC), when asymptotic properties of one-dimensional 
(ID) quantum systems are studied in the thermodynamic 
limit. This is partially because boundary energy correc- 
tions exist under OBC, where eigenstates are not trans- 
lational invariant. Systems with PBC normally contain 
smaller finite size effects, and this property of PBC is 
appropriate for accurate determination of bulk proper- 
ties by means of the finite size scaling. [H-Q 

In numerical studies of lattice models, OBC are of- 
ten chosen for technical reasons. In particular, majority 
of the practical numerical analyses by the density ma- 
trix renormalization group (DMRG) method (^-Q are 
performed under OBC. Concerning to finite size systems 
with PBC, the crucial point in DMRG is the ring-shaped 
geometry, which reduces the decay rate of the singular 
values. Although recent progress in DMRG and the 
tensor product formalisms made it possible to include the 
PBC in a natural manner, numerical implemen- 

tation requires additional computational resources com- 
pared with conventional DMRG analyses. 

A way of suppressing the boundary energy corrections 
induced by OBC is to introduce smooth boundary con- 
ditions (SBC). [O, [13 Recently we proposed a variant 
of the smooth boundary conditions, where the local en- 
ergy scale of TV-site systems is proportional to the de- 
formation function [sinfjTr/TSQ]^ specified by the site in- 
dex 1 < j < - 1. [13|, llj This sine-squared defor- 
mation (SSD) completely suppresses the boundary 
effects when the ground-state energy of a free-Fermions 
model on the ID lattice is considered. In this article we 
generalize the deformation function, which is given by 
(J- = [sin(j7r/iV)]"', where m is a positive integer, [ig} In 



the next section, we examine the effect of this sinusoidal 
deformation by up to m = 5 when it is applied to the 
ID free-Fermion model. It is shown that the case m — 2, 
the SSD, is the most efficient for the suppression of the 
boundary effect. 

Another trial in this article is the application of SSD to 
correlated systems. As typical examples of correlated sys- 
tems, we choose the spinless-Fermion model with nearest- 
neighbor interaction and the extended Hubbard model; 
we report the numerical result obtained by DMRG in 
Sec. III. When the interaction is present, the determina- 
tion of the chemical potential is non trivial. We present 
a systematic way of solving this problem in Sec. IV. We 
summarize results in the last section. 



II. SINUSOIDAL DEFORMATION 

We start from the sinusoidal deformation applied to 
the free-Fermion model on the ID lattice. Consider a 
tight-binding model represented by the Hamiltonian 
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where N is the system size, and t the hopping energy. 
Operators Cj and Cj represent creation and annihilation 
of Fermions, respectively. The parameter a specifies the 
boundary condition, where OBC and PBC correspond 
to a = and a = \, respectively. (The choice a — 
— 1 is known as the anti-periodic boundary conditions, 
which we do not treat in this article.) For each boundary 
condition, the one-particle energy is expressed as 
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where the energy index £ runs from 1 to A^. The ground- 
state energy £'q^'' at half fllhng is obtained as the sum 
of Sf below the Fermi energy e-p — 0. 

Throughout this article we focus on the system size 

AN) 



dependence on the energy per site Eq /N or the energy 

per bond, which is E^J^'^ /N under PBC and is E^q^'^/{N- 
1) under OBC. After a short algebra, one obtains 
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with OBC, where the leading order of the finite size cor- 
rection is proportional to N^^. The correction of the 
same order also exists for the energy per bond. With 
PBC, one finds 
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where the leading correction is of the order of N^'^. The 
difference between Eq. (3) and Eq. (4) chiefly comes from 
the presence of the boundary energy which exists only 
when OBC are imposed. 

The sinusoidal deformation introduces a position de- 
pendent energy scale gj = [sin(j7r/iV)] to each bond 
of the system with OBC, where m is the positive inte- 
ger, [l^ [13] Deforming H'^^'' in Eq. (1), we obtain the 
corresponding free-Fermionic Hamiltonian 
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We have not obtained analytical solution for the one- 
particle spectrum of H^j^^ so far, except for the zero- 
energy state. Thus we perform numerical analyses in the 
following investigations on the ground state. 

Since we are interested in the ground-state energy per 
site (or per bond), we introduce the normalization factor 
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which is the sum of the deformation factors over the en- 
tire system. When m is an odd positive integer, we have 
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FIG. 1: (Color online) Convergence of ep^^ in Eq. (9) with 
respect to A'' at half filling. Data with OBC, PBC, and the 
deformed cases with m = 1 to 5 are shown. 



and when m is an even positive integer, we have 
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We represent the ground-state energy of at half 

filling by the notation i?Q^'* . It is expected that the nor- 
malized energy 
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converges to — 2t/7r in the large N limit in analogy to 
Eqs. (3) and (4). We refer e|,^^ in Eq. (9) to as the 
energy per bond in the following. As a convention, we set 
B(^) = iV - 1 for the system with OBC, and B'-^^ = N 
with PBC, where these values just represent the number 
of bonds. Using this extended definition of B'-^\ we can 
represent energy per bond by Eq. (9) regardless of the 
boundary condition or the presence of deformation. 

We regard t as the unit of the energy in the numeri- 
cal analyses. Figure 1 shows the A^-dependence of gg^' 
in Eq. (9) for the undeformed systems with OBC, PBC, 
and the deformed systems from m = 1 to 5. When the 
PBC are imposed, the convergence of e^f^^ with respect 
to N^'^ is linear, and there is even-odd oscillation with 
respect to the particle number N/2. Similarly the linear 
A'^~^-dependence is observed when m > 2 under the sinu- 
soidal deformation. In the case m = 1, there is additional 
logarithmic correction as shown later. It should be noted 
that when the particle number N/2 is odd, Bq^'' obtained 
with the sinusoidal deformation for m = 2 coincides with 
eQ^"* obtained with PBC. [17| This complete agreement is 
checked down to the smallest digit in numerical precision. 
Throughout this paper we use the exact diagonalization 
in order to reduce any numerical errors to minimum. 
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FIG. 2; (Color online) Asymptotic behavior of ej,'^^ for the relation function (c]c 



FIG. 3: (Color online) Expectation value of the bond cor- 



deformed chains when 1 < m < 5 with respect to eg^^ of the 
system with PBC. A logarithmic correction is present when 
m = 1. 



In order to confirm the N ^-dependence of Cg^'' with 
the sinusoidal deformation under m > 2, we plot the 
difference between eQ^"* obtained with PBC (when N/2 is 
even) and Cg^^ with the sinusoidal deformation. To avoid 

any confusion, let i^pBc ^^^^ E^(2} denote the ground- 
state energy obtained with PBC and with the sinusoidal 
deformation, respectively. We also use a similar notation 



sinusoidal deformation. 



.j+i + cj4.iCj) with respect to j under 



for the normalization factors B 
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TV and b1^2 for the 



normalization factor defined in Eq. (6). Figure 2 depicts 
the magnified difference 
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when N is even. It is shown that the logarithmic correc- 
tion (TVlogA^)"^ is present when m = 1, and is absent 
when m > 2. 

Figure 3 shows the spatial distribution of the bond cor- 
relation function (cjcj_|_]^ -|- c]_|_]^Cj) at half filling when 
N — 1000. The Friedel oscillations induced by the 
boundary are clearly observed when OBC are imposed 
(the asterisks) , and weaker oscillations are observed with 
the sinusoidal deformation when m = 1. Only when 
TO = 2, there are no oscillations at all; we checked the 
uniformity (the translation invariance) of the bond cor- 
relation function down to the 16-digits in numerical preci- 
sion. When m > 3, the boundary effects appear again. In 
this case the bond correlation function toward the system 
boundary does not oscillate, and decreases in monotonic 
manner. Such behaviors for each m might be related to 



the suppression of the boundary corrections in Cq 
Sciency 

SBC proposed in Ref. [lll[l3|- Figure 4 shows the bond 
correlation function for both cases at half-filling, where 
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FIG. 4: (Color online) Comparison of the expectation value of 
the bond correlation function at half filling under SSD (m = 
2) and SBC in Ref. [O, The number M in the case 

of SBC specifies the length of area where the interactions are 
modified near the system boundary. The bottom graph shows 
the numerical details when M = 30 (SBC) with respect to 
SSD. 



the length of boundary area in SBC is chosen as M = 10 
and M = 30 when the system size is iV = 1000. Although 
bulk property is well captured by SBC already for M = 
30, boundary fluctuations are still present. On the other 
hand, the bond correlation function is almost uniform 
away of the boundary. 

Now we discuss the way of treating the deformed sys- 
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tern away from half filling. For the undeformed systems 
with OBC or PBC, it is sufficient to include the chemical 
potential term — mX^^Li into Eq. (f), where n = cjc 



is the number operator. The value of /i adjusts the Fermi 
energy to zero, and is given by 



^i{f) = ~2tCOs{7Tf)., 

where / is the filling factor 
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A natural way of introducing under the sinusoidal 
deformation is to write down the Hamiltonian as a sum 
of the local terms 
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In order to confirm the validity of these constructions 
in Eqs. (11)-(14), we carried out numerical calculations 
for the selected fillings / = 1/4 and / = 1/8. Figure 5 

shows the A^^^-dependence of Gq^'' = E^^^ 



Eq^^ is the ground-state energy for each filling. We plot 
the data only when the particle number p = fN is even. 
Analogous to half filling, the bond energy Cq^'' with PBC 
coincides with that obtained with SSD (m = 2) when p 
is odd. [17] The logarithmic corrections are again present 
when m = 1. 

The uniformity of the ground state obtained with the 
sinusoidal deformation is checked by calculating the oc- 
cupancy (rij) = (c]s) / ^ 1/^' 1/^' ^'^"^ 1/^- Figui'^s 6 
shows (rij) when m = 1, 2, and 3. At half filling, / = 1/2, 
{rij) is always equal to 1/2 by the particle- hole symme- 
try. Even away from the half filling, this uniformity is 
kept when m = 2. There are small fluctuations near the 
system boundary when m = 1 and m — 3. In this way, 
the construction of Hamiltonian in Eqs. (13) and (14) 
away from half-filling is justified, especially for SSD. 
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FIG. 5: (Color online) Convergence of Bq^' with respect to 
at 1/4-filling and 1/8-filling. 
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III. CORRELATED SYSTEMS 

It is expected that the SSD reduces boundary effect 
even when interactions are present between particles. In 
order to check this conjecture, we apply the sinusoidal de- 
formation to correlated systems, by means of the Hamil- 
tonian Eq. (13) written in the linear combination of bond 
operators. We study the system size dependence of the 
ground-state energy Eq^"^ and the uniformity of the sys- 
tem. 



FIG. 6: (Color online) Position dependence of the occupation 
number (rij) = (cjc^) at / = 1/2, 1/4, and 1/8 when m = 1, 
2, and 3 (upper graph). The left and the right lower graphs 
display the occupancy at / = 1/4 close to the boundaries. 



A. Spinless Fermions 

As an example, let us consider spinless Fermions on the 
ID lattice, whose behavior is described by the uniform 
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FIG. 7: (Color online) Energy per bond e^{^^ of the spinless 
Fermions model at V = 1 up to A'' = 16. 
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with PBC, where the system contains the repulsive 
Coulomb interaction V > between neighboring sites, 
in addition to the hopping amplitude t. In this section 
we restrict ourselves to the half-filled case only, therefore 
the chemical potential /i is zero, [l^ The construction of 
the Hamiltonian with OBC and its sinusoidal deforma- 
tion is analogous to what we have done in the previous 
section. The deformed Hamiltonian can be defined by 
putting the bond operator as 



h 



(16) 



and substituting it to Eq. (14). We calculate the ground 
state and the corresponding energy E^^^ at half-filling 
up to = 16. 

Figure 7 shows Gq^^ with respect to N~^. When the 
sinusoidal deformation is applied, we observe the same 
N^'^ dependence as was seen for the non-interacting case. 
In particular when N/2 is odd, Cq^-* obtained for m = 2 
shows a good agreement with the result from PBC. The 
coincidence, however, becomes less accurate with increas- 
ing and there is a deviation around 0.2% in Cq^^ when 
V = 5. Figure 8 shows the occupation number (n^) and 

the bond correlation {cjCj_^_^ + Cj_^_^Cj) when F = 1 at 
half filling. A clear uniformity is observed when m = 2 
as shown in the non-interacting cases. 



FIG. 8: (Color online) Occupation number (rij) and the bond 
energy (ctc^^j -|- ct^j^c^) when V = 1 at half filling. 



B. Extended Hubbard Model 

Now we consider the extended Hubbard model, which 
contains the on-site Coulomb interaction U and the 
neighboring interaction V. In this case there is spin de- 
gree of freedom (cr J,), therefore creation and annihi- 
lation operators, respectively, are represented as cj^. and 
C g.. The bond operator of the extended Hubbard model 



is represented as 



+ (^+lt-5)(^+u-I)] (17) 

To avoid the complexity of determining the chemical po- 
tential fi, we consider the half-filled case where ^ — 
is guaranteed by the particle hole symmetry. Figure 9 
shows the iV~^-dependence of Cg^' for various combina- 
tions of U and V. The coincidence between PBC and 
SSD at half filling occurs when the total number of both 
up- and downs-spin Fermions are odd. 

The occupancy (n^) = {n^^ + n-^) and the bond cor- 
relation function Yj^^'^la'^j+ia + ^1 I'^lf filling 
are shown in Fig. 10. The occupancy (n^) with OBC is 
clearly influenced by the system boundaries, whereas the 
sinusoidal deformations (to > 1) lead to almost constant 
distribution. The case with to = 2 (asterisks) realizes 
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FIG. 9: (Color online) System size dependence of Cg^' of the 
extended Hubbard model at half filling. The full and the 
dashed lines connect the energies Cq'^' obtained by PBC and 
SSD, respectively. 



the minimal position dependence. The bond correlation 
function is influenced by the boundaries in all cases, and 
the position dependence is the weakest when m = 2. 



IV. CHEMICAL POTENTIAL 

So far we have not discussed the proper value of the 
chemical potential away from the half filling when 
the interaction is present. The chemical potential term 

■~f ('^it + ^ii + "j+it + "-j+iJ.) included to 

the bond operator in Eq. ([15]). In the Hartree-Fock 
or the Fermi liquid picture, /i is adjusted so that the 

coincides with 
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particle number p = 

the number of negative-energy quasiparticle states. The 
number of particles is independently represented by the 
derivative of E^^^ with respect to [i. Thus, the relation 
y = —dEQ^^/dfi must be satisfied for the targeted par- 
ticle number p. Within the sinusoidal deformation, this 
relation is slightly modified as 
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according to the position dependence in the energy scale. 

We plot the absolute value \p ~ {—dNe'^^ /d^)\ with 
respect to ^ for the extended Hubbard model in Fig. 11, 
under the SSD. The numerical analysis by exact diago- 
nalization gives /i = —4.1425 for the case iV = 12, p = 4, 
U = 2, and I^ = 1. Figure 12 shows the corresponding 
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FIG. 10: (Color online) Occupation {nj^) = inj{) and the 
bond energy '^„{c]a'^j+-i+c'jj^-^„Cj^) of the extended Hubbard 
model at half filling. 
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FIG. 11: (Color online) The dependence of \p + dNe[^'' /d^Ji\ 
for the two cases / = p/2N = 1/4 and 1/6. 



occupation and the bond correlation functions. There is 
a slight position dependence, since we are dealing with 
relatively small system size with a few particles. The 
position dependence becomes conspicuous if we choose 
either an inaccurate value of /i or we consider the defor- 
mations with m ^ 2. 



V. CONCLUSION AND DISCUSSION 

We have shown that the sinusoidal deformation ap- 
plied to ID quantum Hamiltonians improves convergence 
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ward the thermodynamic hmit, compared to the uniform 
systems with open boundary conditions. Such suppres- 
sion of the boundary effects is confirmed also for interact- 
ing systems, typically for the extended Hubbard model. 

We have not determined the 'canonical' form of the si- 
nusoidal deformation in case where there are long-range 
interactions. It is reported that a small but finite resid- 
ual boundary effect appears in spin chains which includes 
the next-nearest-neighbor interaction, [l^ Application of 
the sinusoidal deformation to higher-dimensional quan- 
tum systems could be a future problem. 

A theoretical puzzle of the sinusoidal deformation re- 
mains in coincidence of the ground-state energy calcu- 
lated with PBC and SSD. The agreement is almost per- 
fect, which strongly suggests a hidden algebraic relation 
between these two cases. 
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